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The nine milhsecond pulsars (MSPs) that have now been detected by Fermi-hNT are providing 
an excellent opportunity to probe the emission geometry of these ancient compact objects. As they 
are radio-loud, one may use the relative phase lags across wavebands to obtain constraints on the 
orientation, size, and location of their radio and gamma-ray beams. We model the gamma-ray 
light curves using geometric outer gap (OG) and two-pole caustic (TPC) models, in addition to a 
pair-starved polar cap (PSPC) model which incorporates the full General Relativistic E-field. We 
find that most MSP light curves are fit by OG and TPC models, while PSPC is more appropriate 
for two others. The light curves of the newest discovery, PSR J0034— 0534, are best modeled using 
outer magnetosphere OG / TPC models of limited extension for both radio and gamma-ray beams. 
We model the radio emission of the other eight MSPs using a fixed-altitude conal model at lower 
altitude. We lastly deduce values for inclination and observer angles (a and C,), as well as the flux 
correction factor, in each case. 
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I. INTRODUCTION 

Our pursuit of understanding gamma-ray pulsars 
has been dramatically furthered by the detection of 
46 pulsars after only six months of observation by 
the Fermi Large Area Telescope (LAT) . Owing to 
the Fermi LAT's superior capabilities compared to its 
predecessor, the Energetic Gamma Ray Experiment 
Telescope (EGRET), a subpopulation of gamma- 
ray millisecond pulsars (MSPs) is now emerging 
twenty seven years after the discovery of the first 
MSP [6^. This population includes PSR J0218-h4232 
(first marginally detected by EGRET (Hi), the first 
i^ermz-detected MSP PSR J0030-h045f, and now a 
ninth MSP PSR J0034-0534 which has nearly phase- 
aligned radio and gamma-ray light curves [3|. It is 
expected that this number will soon grow as Fermi 
accumulates data. Additionally, follow-up radio ob- 
servations of unidentified Fermi sources are also un- 
covering new MSPs [3l| . 

The availability of quality gamma-ray and radio 
MSP data now affords the opportunity to investigate 
the mechanisms and locations of emission in MSP 
magnetospheres. Several different approaches have 
been followed to describe the high-energy (HE) radia- 
tion. Polar cap (PC) models ]Xm, assume near-surface 
emission along with magnetic pair production. This 
establishes a low-altitude pair formation front which 
screens the accelerating electric field [13] ■ Slot gap 
(SG) models 0, ^ envision the formation of narrow 
acceleration gaps in a two-pole caustic (TPC) geome- 
try [l3| . extending from the neutron star (NS) surface 
to near the light cylinder, along the last open B-field 
lines. Lastly, outer gap (OG) models assume that 
photon-photon pair production-induced cascades pro- 
duce HE emission along the last open field lines above 
the null-charge surfaces (where the Goldreich- Julian 



charge density changes sign). Modern extensions to 
this model include a 3D solution [2^ where a small 
positive acceleration E-field reaches to the NS surface. 

We model these ancient, recycled MSPs [8,] using 
3D emission modeling, including Special Relativistic 
effects of aberration and time-of- flight delays, and ro- 
tational sweepback of B-field lines [37J. In addition 
to using geometric OG and TPC models, we find that 
some light curves exclusively favor pair-starved polar 
cap (PSPC) pulsar models, where the pair multiplicity 
is not high enough to screen the accelerating E-field 
and charges are continually accelerated up to high al- 
titudes over the full open-field-line region [l^, [2l|, [sl] . 
For this case, we implement a solution of the accel- 
eration E-field valid up to high altitudes [l^l- Our 
calculation of light curves and flux correction factors 
for the case of MSPs is therefore complementary to an- 
other work j38i] which focuses on younger pulsars, al- 
though our OG and TPC models have non-zero emis- 
sion widths. We describe the details of the various 
models in Section |TT1 and present our results and con- 
clusions in Sections IIIII and HVl 



II. MODEL DESCRIPTION 

Here we briefly describe the models used to gener- 
ate the radio and gamma-ray light curves. A detailed 
description may be found elsewhere j37j . In addition, 
PSR J0034— 0534 presents a special case, and its mod- 
eling will be fully discussed in another publication [36j . 



A. Magnetic Field 

We use an implementatio n of the retarded vac- 
uum dipole B-field solution [l2|, assuming that it rea- 
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sonably approximates the MSP magnetospheric struc- 
ture. This B-field geometry is characterized by an 
asymmetrical PC shape due to rotational sweepback 
of the field lines. We assume constant emissivity [l^l 
along the B-lines in the gap regions of the geometric 
OG and TPC models, while we calculate this explic- 
itly in the PSPC model using the General Relativistic 
E-ficld (Section urn. 

To calculate the direction of photon propaga- 
tion [3^, we find the tangent to the local B-line in 
the co-rotating frame (using the fact that the vacuum 
B-field is approximately the same in this frame as in 
the inertial observer frame), and then aberrate this 
direction via a Lorentz transformation (transforming 
from the instantaneously co-moving frame to the ob- 
server frame) . We also include phase shifts due to ro- 
tational sweepback of the B-ficld lines and travel time 
delays due to the finite speed of light. Our calcula- 
tions are consistent to first order in r jKhc-, with r the 
radial position, and i?LC = c/r2 the light cylinder ra- 
dius, with il the rotational velocity of the MSP. Future 
work may include higher-order corrections, especially 
when considering the force-free magnetospheric struc- 
ture We furthermore explicitly use the curvature 
radius of the B-field lines as calculated in the inertial 
observer frame when calculating curvature radiation 
(CR) losses for the PSPC case. We calculated light 
curves from OG and TPC models for several differ- 
ent gap widths w, while we used w —\ for the PSPC 
model (Section 



B. Accelerating Electric Field 

In the case of the PSPC model, we only consider 
CR losses suffered by electron primaries as they move 
along the B-field lines when modeling the HE emis- 
sion. Previous studies [l^, [2l|, [sH] have used a so- 
lution ^] of the E-field parallel to the B-field {E\\), 
which is valid up to a significant fraction of i?LC7 for 
the PSPC E-field. We now incorporate the small- 
angle-approximation solution of E\ | valid for altitudes 
very close to i?LC l3Q|] • Our implementation of E\ | con- 
serves energy along the relativistic electron primaries' 
trajectories, with electric potential energy being con- 
verted into gamma-ray radiation and particle kinetic 
energy. 



C. Radio Emission Model 

For eight of the Fermi MSPs (excluding 
PSR J0034— 0534), we model the radio emission 
beam using an empirical cone model which assumes 
that the radio emission originates in a cone of 
fixed altitude [111, centered on the magnetic axis. 
We adopt a description [l9| based on fits @ of 
average-pulse profiles of a small collection of pulsars 



at 400 MHz to a core and single cone beam model. 
The emission occurs increasingly close to i?LC as 
the period P decreases (for more or less constant 
P), and is typically located at 10% - 20% of i?Lc 
for these MSPs. Depending on where the observer's 
linc-of-sight intersects the radio cone, the radio profile 
may exhibit zero, one or two peaks. 

In contrast, the near phase alignment of the ra- 
dio and gamma-ray light curves of PSR J0034— 0534 
forces us to depart from the conal framework when 
modeling this newest MSP member. The data seem to 
require caustic origin of the radio (see Section IIII C|) . 



D. Flux Correction Factor 

In order to convert the observed (phase-averaged) 
energy flux Gobs to the all-sky luminosity L-y, which 
defines the gamma-ray conversion efficiency r\^ = 
L^/E-cot, with E-cot the spin-down luminosity, we need 
to calculate a flux correction factor (/n). This fac- 
tor accounts for the fact that an observer only sees a 
small part of the total radiation: that coming from a 
slice through the emission beam. The flux correction 
factor is defined as follows 1381: 



obs ; 



JJ Fj{a, (j)) sin QdCdcj) 
2/F^(a,CB0)d0 ' 



(1) 



(2) 



with F^ the photon flux per solid angle ('intensity'), 
and C-E the Earth line-of-sight (assuming similar dis- 
tributions of gamma-ray photon and energy fluxes in 
(C,0)-space, i.e., F^{Ci4>)/F^, tot ~ G^{(,(l))/G^^tot, 
with <j) the phase). Results of calculations of /ji for 
different pulsar models appear in Table [J 



III. RESULTS 

We generated light curves for each of the different 
pulsar models, using a range of inclination and ob- 
server angles (a = C = 5° — 90°, in 5° intervals), pe- 
riods P = 2, 3, and 5 ms, and gap widths w (ranging 
from 0.05 - 0.4 for the OG / TPC models, and w = 1 
for the PSPC case). We matched model light curves 
to the MSP gamma-ray and radio data by eye (al- 
though statistical uncertainties in the data as well as 
similar predicted light curves shapes for similar ranges 
of a and C may complicate unique matching). Fig- 
ure [T] shows an example of this for four MSPs. We 
did not find any satisfactory fits for the geometric PC 
models or for OG and TPC models with w = 0. In 
addition, the radio cone model fits the data (exclud- 
ing PSR J0034— 0534) quite well overall, except for 
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TABLE I: Inferred values of a, (, and fn{a,C,P) 



PSR Name 


otpc 


Ctpc 


QOG 


Cog 


apspc 


Cpspc 


Q^radio 


Cradio 


Ref. 


/n.TPC 


/n,OG /n,pspc 


J0030+0451 


70° 


80° 


80° 


70° 


— 


— 


~ 62° 


~ 72° 


[26] 


1.04 


0.90 


— 


J0034-0534 


30°" 


70°" 


30°' 


70°' 


— 


— 


— 


— 




0.74" 


0.45* 


— 


J0218+4232 


60° 


60° 


50° 


70° 


— 


— 


-8° 


~ 90° 


[34J 


1.06 


0.63 




J0437-4715 


30° 


60° 


30° 


60° 






20° - 35° 


16° - 20° 


[17, 28] 


1.23 


1.82 




J0613-0200 


30° 


60° 


30° 


60° 






small f3 




[39] 


1.19 


1.76 




J0751+1807 


50° 


50° 


50° 


50° 












0.80 


0.65 




J1614-2230 


40° 


80° 


40° 


80° 












0.83 


0.64 




J1744-1134 










50° 


80° 












1.19 


J2124-3358 










40° 


80° 


20° -60° (48°) 


27° - 80° (67°) 


[27] 






1.29 



"Limited TPC model. 
''Limited OG model. 



the case of PSR J0218+4232, which seem to require a 
wider cone beam. The MSP radio beam may be rela- 
tively large, as its total size and annular width scales 
as p~*'-35 in our model. 

From our light curves fits, we infer values for a and 
C for each MSP. These are summarized in Table [I] 
(labeled with subscripts 'TPC, 'OG', and 'PSPC'), 
and compared with values obtained from radio polari- 
metric measurements (labeled with subscripts 'radio'). 
Our calculated values for /n(a, C, P) indicate that this 
factor is typically of order unity for the best-fit geome- 
tries we consider here, with the OG model usually 
predicting lower values than the TPC model. Our re- 
sults qualitatively resemble other calculations for OG 
/ TPC models applied to young pulsars [11], while 
the difference encountered when using the unscreened 
PSPC model (vs. the PC model) reflects the funda- 
mental physical dissimilarity of the magnetospheric 
structure of pair-starved MSPs and younger pulsars. 



A. OG / TPC Case 

Both OG and TPC models have a preponderance of 
double-peaked light curves at similar phases. Sharp, 
solitary peaks exist for some regions in phase space in 
OG models, while the corresponding TPC-peaks usu- 
ally have additional low-level features. No emission 
originates below the null charge surface (C = 90°) 
in the OG model, implying that emission is only 
visible from one magnetic pole, in contrast to the 
TPC model where emission is visible from both poles. 
Consequently, OG models do not exist at all angle 
combinations, while TPC models do. (We impose a 
maximum emission radius rmax — 1-2-Rlc for these 
models, together with a maximum cylindrical radius 
Pmax = 0.95i?LC, while we assume that the emission 
starts at the stellar surface r = R.) 

Six MSPs (all in Table |T] excluding 
PSR J0034-0534, PSR J1744-1134, and PSR 



J2124-3358) are well fit by OG and TPC models. 
For these cases, the gamma-ray light curve lags the ra- 
dio. In addition, PSR J0030-h0451, PSR J0218-h4232, 
and PSR J 1614- 2230 (as well as PSR J0034-0534; 
see Section IIIIC|) exhibit double-peaked light curves, 
which signifies the existence of screening electron- 
positron pairs that form the OG or TPC emitting 
structure. 



B. PSPC Case 

For the PSPC model, we find mostly single-peaked 
gamma-ray profiles roughly in phase with the radio 
(for single radio peaks), which increase in width when 
P decreases, especially for the radio. For large im- 
pact angles (3 — C, — a, only gamma-ray radiation is 
visible, since the gamma-ray beams are believed to 
be larger than the radio ones (possibly explaining the 
phenomenon of 'radio-quiet' pulsars 2]). 

Double-peaked profiles occur only for large C, (not 
for both large a and C as in a constant-emissivity PC 
model), because the presence of 'favorably-curved B- 
field lines' at large a is reflected in the azimuthal de- 
pendence of Emission is therefore usually only 
visible from a single magnetic pole, as in the OG 
case. The gamma-ray emission originates on all open 
field lines, while the radio originates at a fixed height. 
The gamma-ray profile therefore appears at earlier 
phases than the radio (with the caustic peaks being 
washed out), so that the radio peak lags the gamma- 
ray peak. This is the case for PSR J1744-1134 and 
PSR J2124-3358, which are exclusively fit by the 
PSPC model. For these MSPs, the gamma and ra- 
dio emission originate from the same magnetic pole, 
well above the NS surface. 
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FIG. 1: Plots of the observed and fitted gamma-ray and radio ligiit curves for several MSPs: PSR J0030-I-0451, PSR 
J0218-I-4232, PSR J0437-4715, and PSR J1744-1134. In each case, in the upper panels the histograms represent the 
Fermj-LAT data above 0.1 GeV [3, the horizontal dashed line the estimated background level, the green lines are 
OG fits, and magenta lines are TPC fits (see Table |l]). The bottom panels indicate radio light curves, where blue lines 
represent the radio data, while the magenta and green lines correspond radio model fits with the same (a, C,) combinations 
as those of the respective OG and TPC fits. (For the last two radio panels, the OG and TPC models are for the same a 
and 1^, giving a unique radio conal beam fit.) Light curves were chosen to show different classes of gamma-ray models: 
panel (a) and (c) indicate OG / TPC models with large a, panel (e) shows OG / TPC models with small a, and panel (g) 
shows a PSPC model. 



C. Limited OG / TPC Case 



similar to what is observed in the Crab pulsar. 



We find that a combination of a radio PC model 
and a gamma-ray OG or TPC model does not repro- 
duce PSR J0034-0534's light curves [36i]. Rather, a 
co-located caustic origin of both radio and gamma-ray 
emission is probably required The so-called 'lim- 
ited OG / TPC models (Figure ^ refer to OG and 
TPC models where we have limited the extent of the 
emission regions along the B-field lines, with the radio 
subsumed into the more extended OG / TPC regions. 

The minimum and maximum emission radii are 
quite well constrained by the profile shapes and radio- 
to-gamma phase lags. We found good fits for a ~ 
30°, C = 70°, w = 0.05, and using (r,„i„, Tmax) = 
(0.12, 0.9)i?LC for the gamma-ray profile, while assum- 
ing (ri„in,rmax) = (0.6, 0.8)i?LC for the radio profile. 
This is the first time that near phase alignment of 
the radio and gamma-ray pulses is seen for an MSP, 



IV. DISCUSSION AND CONCLUSIONS 

The gamma-ray observations now support a 
stronger 'familial' link between the MSP and young 
pulsar classes in the sense that we observe spectral and 
temporal characteristics in their HE emission which 
are very much alike. This indicates the operation 
of very similar emission mechanisms regardless of the 
aeons separating the different pulsar generations. In- 
terestingly, it is noted 2] that the MSPs and young 
pulsars observed by Fermi share magnetic fields at 
the light cylinder which are comparable in magni- 
tude, possibly providing further corroboration for this 
idea (see also [lOl)- Also, explanation of the surpris- 
ing inference that there must be a significant num- 
ber of pairs in MSP magnetospheres to support the 
OG / TPC framework, may require a B-field struc- 
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FIG. 2: Data and model light curves for PSR J0034— 0534. The top histogram indicates Fermi-LAT data above 
0.1 GeV 0] along with limited TPC ( green) and limited OG (red) fits. Similar for the bottom panel, where the solid 
blue line is the Westerbork Synthesis Radio Telescope (WSRT) 324 MHz radio profile. We used a = 30°, C = 70°, and 
w = 0.05 for all light curves. 



ture departing significantly from the usual dipole pic- 
ture U H^. Thus, the surface B-fields of younger 
pulsars and MSPs may be of similar magnitude (un- 
less the younger pulsars' B-fields are also boosted, but 
by a different mechanism than in the MSP case) , and 
the low-altitude B-field structure should play an im- 
portant role in the magnetic pair creation physics en- 
countered in TPC and PC models. Additionally, re- 
cent measurements indicate large MSP masses of up to 
^ 1.7 Mq [l^, larger than the average mass of young 
pulsars, and this may impact MSP B-field properties 
and enhance pair creation probability. 

We were able to demonstrate a bifurcation in the 
HE MSP population, distinguishing between MSPs 
with light curves modeled by the OG / TPC models on 
the one hand, and those fit by the PSPC model on the 
other hand. This is due to the fact that we modeled 
both the gamma-ray and radio curves for each pul- 
sar. Both the shapes and the relative radio-to-gamma 
phase lags provided by the data could then be used 
to constrain the best-fit model type and geometry for 
each MSP. 

However, PSR J0034-0534 respresents a third MSP 
subclass, as we inferred that an outer magnetosphere 
caustic origin of both the gamma-ray and radio emis- 
sion is required to explain its nearly phase-aligned 
light curves. We therefore conclude that the HE emis- 



sion must come from the outer magnetosphere in all 
MSP models considered, as is also believed to be true 
for the bulk of the gamma-ray pulsar population Q • 

Our calculations of the flux correction factor imply 
a wide beaming angle, which derives from the fact 
that we obtain best fits for large impact angles. The 
larger radio beam widths of MSPs compared to those 
of canonical pulsars furthermore suggests that there 
should be relatively few radio-quiet MSPs. 

Studies involving phase-resolved spectroscopy of 
Fermi-LAT data, as well as the anticipated discovery 
of several more gamma-ray pulsars, should make for 
a vibrant field of research for the foreseeable future. 
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